RECOVERING THE TWISTING FUNCTION IN A TWISTED WAVEGUIDE 

FROM THE DN MAP 



MOURAD CHOULLIf AND ERIC SOCCORSIf 

Abstract. We consider the inverse problem of determining the twisting function in a infinite cylindrical 
twisted waveguide from the corresponding DN map. This problem, which is naturally linked to some in- 
verse anisotropic conductivity problem in a straight waveguide, remains generally open, unless the twisting 
function is assumed to be affine. Namely we prove Lipschitz stability in the determination of affinc twisting 
functions from the DN map. This result still holds true upon substituting a suitable approximation of the 
DN map, provided the first derivative of the twisting is sufficiently close to some a priori fixed constant. 
Key words : Dirichlct Laplacian, twisted infinite cylindrical waveguide, twisting function, DN map, sta- 
bility estimate. 
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1. Introduction 

Let a; be a bounded domain of M 2 . To 51 = lu x R and 9 G C 1 (R) we associate the infinite twisted 
cylindrical domain 

= {(Re( X3 )x',xs); x' = (xx,x 2 ) G u, x 3 G R}, 
where denotes the rotation in M 2 of angle £ G R. Twisted waveguides modeled by f2g exhibit interesting 
propagation properties such as the occurence of propagating waveguide modes with phase velocities slower 
than those of similar modes in a straight waveguide. This explains why twisted waveguides are at the center 
of the attention of many theoretical and applied physicists (see e.g. [Kal |DRl KEl INZG| [Shi IWi] lYMj ). 
Moreover, it turns out that they are the source of challenging spectral and PDE problems, some of them 
having been extensively studied in the mathematical literature (see e.g. |BK[ lEKK , KK, KSl IKZll IKZ2j ). 

Nevertheless, the inverse problem of identifying the twisting function from the Dirichlet-to-Neumann (DN 
in short) map, has, to our knowledge, not been examined in this framework yet. The study of this open 
problem is actually the main purpose of the present article. Namely, we consider in this paper the following 
boundary value problem for the Laplacian in the twisted waveguide fig, 

Av(y) = 0, y G Qg, 
v(y) = g(y), y e dn e , 

and we address the problem of recovering the twisting function 9 (actually its first derivative) from the DN 
map 

(1.2) l e :g^B{y)Vv{y)- V (y), 

l 



(1.1) 
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where 

A -y 2 9'(y 3 ) 
B(y) =01 yi 9'(y 3 ) 
\0 1 

More specifically, our aim is the stability issue for the problem of determining 9 from Ag. 

The DN map Ag acts on functional spaces depending on 9. Therefore it is not well suited to the analysis 
of this inverse problem. This difficulty may be overcome by turning (|l.ljl - ()1.2|l into an equivalent system 
associated to some ^-independent DN map. To this purpose we introduce 

* - (* £ : 

and put u(x) = v{Tg^ X3 ^{x' ,x 3 )), x = (x',x 3 ) G f2. By performing the change of variable y = Tgr X3 )(x' ,x 3 ) 
in we find by direct calculation that u is the solution to the following boundary value problem with 

an elliptic operator in the divergence form 

n „x f div(A{x',9'(x 3 ))Vu) = 0, x G Q, 

( j \ u(x) = /(a?), x G dSl, 

where, f(x) = g(Tg^ X3 \(x' , £3)), £ G 90, and the matrix ^4 is given by 

(l + x 2 t 2 -x 2 xit 2 -x 2 t\ 

~x 2 xit 2 l + xft 2 xit , x'Gw, set. 

— a^i £ii 1 / 

Moreover, it holds true that 

(1.4) S( 2 /)Vu(y)-^(2/)=^(x / ,^(a;3))Vu(a;)-^(a;), y = r fl(a ,, ) (s / ,« 3 ) 1 a: = (x', x 3 ) G fi. 

This identity justifies the choice of the boundary operator B appearing in ()1.2j) . which provides the appropri- 
ate Neumann condition, given by the right hand side of (|1.4p . on u. Formally, p.4[) indicates that recovering 
9 from A# is the same as determining 9 from the following DN map 

A e : / ->• AVu ■ v. 

The major part of our work will therefore be devoted to studying Ag in view of establishing stability in the 
identification of 9 from Ag. In light of (|1.3|l - (|1.4|) we notice that this is the same kind of inverse anisotropic 
conductivity problem, but stated here in an unbounded domain, as the one studied in a bounded domain 
by Alessandrini [X] and Alessandrini and Gaburro jAGlj . |AG2j (see also Gaburro and Lionheart |GL| ). 
However, it turns out that the usual monotonicity assumption on the conductivity, which is essential to the 
identification of A from the DN map in this approach, is not fulfilled by the matrix A under consideration. 
This is the main reason why the inverse problem associated to (|1.3p remains open for general twisting 
functions 9 G C 1 (M). Nevertheless, in this paper we are able to prove Lipschitz stability in the determination 
of 9 from the DN map when 9 is cither an affine function or sufficiently close to an arbitrarily fixed constant. 

The paper is organized as follows. Section 2 contains the definition and the main properties of Ag, needed 
in the proofs of the coming sections. The first part of section 3 explains why the approach developpcd 
by Alessandrini and Gaburro in [AGlj and |AG2j does not apply to the problem under consideration. 
Nevertheless, we prove in the second part of section 3 that twisting functions close enough to some arbitrarily 
fixed constant may actually be identified by solving an inverse conductivity problem which satisfies the above 
mentioned monotonicity condition. Further, the particular case of affine twisting functions is examined in 
section 4. This is by means of the Fourier transform with respect to the variable x 3 , in order to bring the 
original problem into some anisotropic conductivity problem stated in oj. The corresponding conductivity 
matrix satisfies a weak monotonicity condition, allowing to claim stability in the determination of the twisting 
function from some suitable DN map. Finally, the original DN map Ag is rigorously defined and linked to 
Ag in section 5. 
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2. The DN map 

Extension lemma. As we are dealing with an infinitely extended domain f2, we start by proving the 
following useful extension lemma. 

Lemma 2.1. Let g G £P +1 / 2 (R; H s (du>)) for s = 3/2 or s = 1/2. Then there exists G £ £P +1 / 2 (R; H s+1 I 2 {uj)) 
such that G(t) = g(t) on duj and 

(2-1) l|G|| fl - s +i/2( R;fl - s +i/2( tJ )) < C(u})\\g\\ H s+-L/2(p l . H ,( L d u fl, 

where C{uj) is a constant depending only on lu. 

Proof. We detail the proof for s = 3/2, the case of s = 1/2 being treated in a similar way. Let us first assume 
that g G Q°(R; H 3 / 2 {duj)). For each h 6 H 3 / 2 {dw), the boundary value problem 



(2.2) 



AH = in u, 
H = h on du, 



admits a unique solution H G H 2 (uj), according to |LM) . Moreover there exists a constant C(u>), depending 
only on u, such that the following estimate holds true: 

(2-3) \\H\\ m(u) <C(Lj)\\h\\ HV2idu>) . 

Let G(t), t G R, be the solution of (|2.2[) corresponding to h = g(t). Since G(t) — G(s), s, t G R, is the 
solution to the boundary value problem (|2.2p with ft, = g(i) — g(s), we deduce from (|2 . 3[) that 

\\G(t) - G(s)\\ H2(u) < C(uj)\\g(t) - g(s)\\ H3/2(du) . 

Therefore G G C(R; H 2 (u)), and we have 

(2-4) ||G|| L 2( R . // 2( tlj)) < C(Ct>)||ff||z,2( K ;ff3/2(8u>))! 

from (|2.3[) . Let us next consider the solution A'(t) to the boundary value problem (|2.2[) associated to 
/i = g'(t). Similarly we have K G C(R; H 2 (uj)), and (gjj yields 

||G(t + «) -G(i) - S A(t)|| H 2 M < C(w)\\g(t + s) - g(t) - sg'(t)\\ H3/Hdu)) , 

since G(i + s) — G(i) — sK(t), t, s G R, is the solution of (|2.2|) associated to /i = g(t + s) — g(t) — sg'(t). 
From this then follows that G G C 1 (R] H 2 (uj)), G' = K, and 

(2-5) ||G ! '|| i 2( R . // 2( tlj)) < C{uj)\\g'\\ L 2 { f L . H 3/2 {dbj) y 

By substituting G' for G in the above reasoning, we get that G G G 2 (R; H 2 (oj)), and 

(2-6) l|G"|| L 2( R . // 2( w)) < C(u;)|[flr // || i a( R . J y3/2( aw )). 

Putting (|2.4 |) -([2~5 ]1 together, we thus find out that 

(2-7) IIGIIiT^Rjif 2 ^)) < C( w )ll5ll_f/ 2 (K;ff 3 / 2 (aw))- 

Further, g G H 2 (R; H 3 ^ 2 (8oj)) being fixed, we consider a sequence (f/„) n in C^°(R; H 3 / 2 (du))) converging 
to g in iJ 2 (R;ff 3 / 2 (<9w)). By uniqueness, G n — G m , where G„ (resp. G rn ) denotes the corresponding 
extension of g n (resp. <?,„), extends g„ — g m to iJ 2 (R; H 2 {oj)). Therefore, (|2.7[) yields 

|G„ - G m ||H 2 (R;_f/ 2 (Lj)) < G(w)||5„ - gm\\H 2 (M;H 3 / 2 (duj))> 

so (G„)„ is a Cauchy sequence in H 2 (M.; H 2 (uj)) . We call G its limit. From the continuity of the trace 
operator 

W G H 2 {R;H 2 (uj)) -> VF| 9f2 £ H 2 (R; H 3/2 (dco)), 

we see that G extends 5. Thus we end up getting (|2. 1[) by applying (|2.7[) for j„, n £ N, and sending n to 
infinity. □ 
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Fix s = 3/2 or s = 1/2. Since the trace operator w G H s+1 / 2 (uj) — > W\q u G H s {dui) is bounded, then the 
same is true for 

G G H s+1 ' 2 {M.-H s+1 ' 2 {uj)) -+ G\ 9n G F S+1 / 2 (R; H s (duj)). 

Hence, 

\\\9\\\H°+ 1 / 2 (R;H°(duj)) = mi {l|G , ||i/ s +i/2( R . ffs +i/2( tJ )); G = g on <9S1}, 

is a norm on iF +1 / 2 (R; H s (dui)), which is equivalent to the usual one. In the sequel, we shall use either one 
of these two equivalent norms, each of them being denoted by the same above mentioned symbol. 

Solution to the boundary value problem (|1.3|) . Let us first make the following remark on the uniform 
cllipticity of A, where A denotes either A(x',t) or A(x' , #(2:3)), as defined in the previous section. For all 
C G K 3 , x' G ui and t G R, we have 

A(x',t)( • C = Ci 2 + Cf + Cl - 2te 2 CiC 3 + 2teiC 2 C 3 + t 2 (x 2 Ci ~ a;iC 2 ) 2 

= Ci 2 + Cl + (Cs + t(x 2 (i - X1C2)) 2 , x' = (xi,x 2 ) ew,ieM ; 

by a straightforward computation. For every x 1 G u> and t S R, this entails that A(x', t)( ■ ( = if and only 
if C = 0. Therefore, since uJ x [t,t] is compact for all real numbers t < t, there exists A > 1, depending on u>, 
t and t, such that we have 

(2.8) A- 1 |C| 2 < A{x\t)C-C< A|C| 2 for all x' G w, i G (el 3 . 

In order to define the DN map associated to the boundary value problem (|1.3[) , we first need to solve this 
later. To this end, pick / G ff x (R; H 1 / 2 {duj)) and F G H\n) such that F = / on <9ft, where H^R, H 1 ^)) 
is identified with if 1 (SI). Notice that the existence of such a function F is guaranteed by Lemma [3.11 In 
light of (|2.8p and the Lax-Milgram lemma, there is a unique v G -ffg(Sl) solving the variational problem 

(2.9) / AVv ■ Vwdx = - AVF ■ Vwdx, for all w G H£(Q). 

Jn Jn 

Hence u = v + F is the unique weak solution to the boundary value problem (|1.3|) . That is, u satisfies 
the first equation in (|1.3j) in the distributional sense and the second equation in the trace sense. Moreover, 
taking w = v in ()2.9[) . we obtain from Poincare's inequality (which holds true for SI since uj is bounded) that 
IMIi? 1 ^) < Cll-^llff^a) f° r some constant C > depending on u. Therefore we have 

\\u\\ HH n) < C\\F\\ HHa) , 

where C denotes some generic positive constant depending on u>. Finally, as F may be chosen in if 1 (SI) so 
that F = f in dfl, we deduce from this and (|2.ip that 

( 2 -10) IMIifi(n) < C||/lli?i(R;HV2(a w ))- 

Definition of the DN map. Prior to defining the DN map we need some technical result stated in the 
coming proposition. To this purpose we preliminarily introduce the following H (div)-type space, 

F(div A ,0) = {P G £ 2 (S1) 3 ; div(AP) G L 2 (n)}, 

and recall that the dual space of H 1 (R; H 1 / 2 {dui)) is denoted by H~ l {R\ H- 1 / 2 (du})). 

Proposition 2.1. Let P G H(div A , SI). Then AP ■ v G ff-^R; H~ 1 / 2 {du))) and 

(2.11) \\AP ■ u\\ H - 10ltH - 1/a[Bu)) < C(||P|| ia(0 ) + ||div(AP)|| ia( n))0 
In addition, the following identity 

(2.12) (AP-v,g) = I Gdiv(AP)dx+ / AVG ■ Pdx, 

Jn Jn 

holds true for any g G iF(R; H 1 ^ 2 {dw)) and G G iF(Sl) such that G = g on dfl. Here (•,•) denotes the 
duality pairing between H X (R; H 1 / 2 ^)) and its dual F'^R; H-^ 2 (duj)). 

x This means that the operator P 6 Cg°(Q) —>■ AP ■ v £ C°°(dQ) can be extended to a bounded operator from H(divA, 
into H~ 1 (R;H- 1 / 2 (dLu)). 
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Proof. We first consider the case of P E C$°(n) 3 . Fix g G H 1 (R; H 1 / 2 (dixj)) and let G G if 1 (ft) be chosen 
in accordance to Lemma |3. II so that G — g on dCl. Since P has a compact support, we have 

(2.13) / Gdiv(AP)dx = - [ AVG ■ Pdx + [ gAP-vda, 

Jn Jn Jan 



from Green's formula, whence 



gAP ■ vdo 

an 



< C\\G\\ H i {n) (\\P\\ LHn) + ||div(AP)|| i2(a) ) 



By taking the infimum over {G E H 1 ^), G = g on dft} in the right hand side of above estimate, we find 
that ([2TTT) holds true for every P G C^(Jlf. 

Further, pick P G H{div A , fi). The set C£°(n) 3 being dense in H(div A , fi), as can be seen by mimmicking 
the proof of Theorem 2.4 in |GR| . we may find a sequence (Pk)k in C^°(n) 3 converging to P in iJ(divA, fi). 
Moreover, due to p. lip . (AP k ■ v) k is a Cauchy sequence in H~ 1 (R; H~ 1 / 2 {duj)). Therefore (AP k ■ v) k has 
a limit in H~ 1 {R; i7" 1/2 ((9w)), which is denoted by AP ■ v, and (pTH]) follows readily from (f2T3|) . □ 

Let u denote the i? 1 (J7)-solution to (jl.3j) . By applying Proposition 12 . 1 1 to P = Vw, we deduce from (|2.10|) 
that 

A e : / -)• AVw • f 

is well defined as a bounded operator from 7J 1 (M; H x / 2 (duij) into _ff~ 1 (M; H~ 1 ' 2 (dcj)). Moreover the follow- 
ing identity 

(2.14) (A e f,g)= [ A\7u-VGdx, 

Jn 

holds true for all g E H^R; H^ 2 (duj)) and G E H^fl) such that G = g on 9fi. 

Further, by taking G = v in (|2.14j) . where u is the solution to ([1.3)1 with / replaced by g, we find out that 



(Agf,g)= / AVu-Vvdx = / Vu • AVvdx. 
Jn Jn 



Therefore we have 

(2.15) (A e f,g) = (f,A e g), for all f,g G H^H 1 ' 2 ^)). 

This proves that -^e | jji (g j) = wncre H 1 (R; H 1 / 2 (doj)) is identified with a subspace of its bidual 
space. 

Finally, for i = 1,2, put Ai = A(x', 6i(xs)) and Aj = Ag i , and let Ui E 7J 1 (fi), i = 1, 2, be a weak solution 
to the equation 

div(AiVwj) = in CI. 
By applying (|2.14[) to / = u^ dQ and g = W3-i| an , i = 1, 2, we get that 

(AiUx,U2) = / A\\7u\ ■ Vu 2 dx and (A2W2, Mi) = / A2VW2 ■ Vwidcc. 
In light of CHS]), this yields: 

(2.16) ((Ai - A 2 )iii,tta) = I (Ai- A 2 )Vm ■ Vu 2 dx. 

Jn 

Restriction of the DN map. We turn now to establishing some smoothness property for the restriction 
of Ag to H 2 (R; H z / 2 {doj)). We shall assume for this purpose that f2i = lo x (—1,1) has if 2 -regularity 
property, i.e. that for every / E L 2 (Cl) and any matrix-valued function C = (CV, (x))i<i.j<3 with coefficients 
in W 1,co (Cli) satisfying the ellipticity condition 

3a > 0, C(x)£ • £ > a|C| 2 , for all C G C 3 , x G Oi, 
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the following boundary value problem 

div(CVu;) = / in fii, 



w — on dili, 

has a unique solution w G H 2 (fli) obeying 

\\w\\ H 2 {ni) <C(a,M)\\f\\ LHnih 

for some constant C(a, M) > depending only on a, M — maxKijo an d UJ - 

We notice that f2i has iJ 2 -regularity property if and only if fi a = cj x (—a, a) has if 2 -regularity property 

for some a > 0, and from jGrj that fii has -ff 2 -regularity property provided u is convex. 

Having said that we may now prove, upon identifying H 2 (Sl) with H 2 (R; H 2 (uj)) , the following claim, 

which is a cornerstone in the derivation of smoothness properties for the restriction of Ag to H 2 (JS.; H 3 ^ 2 (duj)). 

Theorem 2.1. Assume that 8 G C 1,1 (M) and that f2i has H 2 -regularity property. Then for any f G 
_ff 2 (R; H 3 / 2 (dui)), the boundary value problem (jl.3|) has a unique solution u G H 2 (Q). Moreover if\\0\\ci,ing.) < 
M, /or some M > 0, we may find a constant C > 0, depending only on M and ui, such that we have: 

( 2 - 17 ) II«IIh 2 (0) < C||/llff2(M;i? 2 /3(g a; )). 

Proo/. For M > fixed, let 6» G C 1 ' 1 ^) be such that ||0|| c i.i (R) < M and pick / G H 2 (R; H 3 / 2 (du>)). Wc 
know from Lemma T3. II that there exists F G H 2 (R; H 3 / 2 (ui)) such that F = f on 9f2 and 

(2-18) ll^ll^(n) < C(u)\\f\\ 

Set 

* = div(AVF). 

In light of ()2.8[) there is a unique uq G i?g(Q) satisfying simultaneously 



/ AVu ■ Vvdx = i Wvdx, for all v G Hq(Q), 
Jn Jn 



(2.19) 
and 

(2.20) \\uo\\m(Q) < 0)11*1^2(0), 

for some constant Co > depending on w and M . 

Further, wc consider £„ G Cq°(— (n+ l),n+ 1), for n > 1, such that £„ = 1 in [— n, n] and ||£,'J|oo < 1/2, 
ll^nlloo < 1/2- For every v G /f^fi), we get using standard computations that 

/ AV(£ n w ) • Vvdx = / AVuq ■ V(£ n v)dx - / (AVu ■ V£, n )vdx + / (AV£„ ■ Vv)u dx. 
Jn Jn Jn Jn 

An integration by parts in the last term of the right hand side of the above identity providing 

/ (AV£„ ■ \7v)uodx = — I (j4V£„ • \7uo)vdx — I div(AV£ irl )uovdx, 
Jn Jn Jn 

we next find out that 



AW(£,„u )-Vvdx = / AWu Q -V(£, n v)dx- I (A\7u -\7^ n )vdx~ (A\7^ n -Vu )vdx- div (AV£ n )u vdx. 
Jn Jn Jn Jn 

Since A is symmetric, it follows from the above equality (|2.19[) that 

/ A\7(£ n u ) ■ Vvdx = / y^vdx -2 (AV£„ • Vu )vdx - / div (AV£ n )u vdx, for all v G H^(Q). 
Jn Jn Jn Jn 

As a consequence the function £„uo G i?Q (Q„_|_i), where we recall that Q, a = uj x (—a, a) for any a > 0, is 
the solution to the variational problem 

(2.21) / AV(£„w ) -Vvdx = tyvdx, for all v G H^(Q n+1 ), 
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with 

* = - 2AV£„ • Vw - div(AV£„)«o- 
The next step of the proof is to make the change of variables (x 1 , x 3 ) E — s> (a/, 2/3) = (x', l/(n + l)x3) S 
12 1 in (f2~2Tj) . To this end, we set 

/l \ 
J„ = 1 , n > 1, 
\ 1/n / 

and for all (x',y 3 ) £ f2i, we introduce the following notations: 
A(x', 2/3) = l/(n + 1) J„+iA(x', (n + 1)2/3) Jn+x, 
i(ya) = &»((» +l)lfo), 
u{x',y 3 ) = u (x', (n + 1)2/3), 
w„(x', 2/3) = £ n ((n + l)y 3 )u (x', (n + 1)2/3), 

div(P(x',2/3)) = d Xl Pi(x',y 3 ) + d X2 P 2 (x',y 3 ) + l/(n + l)d y3 P 3 (x' ,y 3 ), 

£(x', 2/3) = l/(n + 1) [tf (a/, (n + l)y 3 ) - 2J n+1 A(a/, (n + l)y 3 ) J„ +1 V^ 3 ) ■ Vu^.i/a) 

- djv(A(a:', (n + l)y 3 ) J n +i V|(j/3 3/3) • 

By performing the above mentioned change of variables, we find out that w n E i?o (^x) is the solution of the 
variational problem 

/ AVw n -\7vdx= / tyvdx, for all v E -ffg (Oi). 

As r^i has _ff 2 -regularity property by assumption, we get by straightforward computations that llifnlliWQi) < 
C(M, w)||*|| L 2 (f2l) . On the other hand, since ||*.]| L 2 (ni) < (n + 1)~ 3 / 2 C(M, w)||*|| L 2 (n) , from (12T2"0T) . it holds 
true that 

(2.22) IKIIflapio < (n + l)- 3 / 2 C(Jkf,a;)[|*|U 2( o). 
Moreover, putting (|2.20[) and (|2.22l) together we obtain that 

(2.23) UnU \\ H z(Q) < C(M,«)||*|| £a( n). 

Therefore, upon eventually extracting a subsequence, we may assume that (£ n Uo)n converges weakly to u in 
H 2 (fl). On the other hand, since (£ ra wo)n converges to uq in L 2 (f2), we have uo E H 2 (tt) and thus (£, n uo) n 
converges weakly to uo in H 2 (fl). Further, the norm || • |_f/2(Q) being lower semi-continuous, we have 

||uo||jya(n) < liminf ||^ || ff 2 (n) < C(M, w)||*|| L 2 (n) , 

from (j2.23|). Bearing in mind that ||*||i,a(n) < C(M,u))\\F\\ H ^(n) and \\F\\ H 2 (n) < C{u)\\f\\ H 2 {R . H 3/2 (du]) , 
this entails that 

||wo||fl-2(n) < C(M,u})\\f\\ H 2(^. H s/2^). 
As a consequence u = uo + F E H 2 {Vt) is the unique solution to (|1.3[) . and it satisfies in addition 

H|jT*(£)) < C( A A^)ll/llH2(K;ff3/2 (aw)) . 

This proves the result. □ 

For each 9 obeying the assumptions of Theorem 12.11 the mapping 

A e : / E H 2 {R; H 3 ^(duj)) -¥ d v u E H 2 (R; H^ 2 (du)), 

where u denotes the unique if 2 (f2)-solution to (|1.3[) . is well defined by Theorem 12. II 

Further, bearing in mind that Cq°(R; H 2 {lo)) is dense in H 2 (il) and arguing in the same way as in the 
derivation of Lemma 13. 1[ we find out that the trace operator 

t : w E H 2 (Q) -> d v w E H 2 (R;H 1/2 (duj)) 
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is bounded. Moreover, we have ||r|| < ||r||, where r : w G H 2 {lu) — >• d u w G iJ 1 / 2 (9w) denotes the usual 
trace operator. From this and (|2.17[) then follows that ||Ae|| < C(M,w), as a linear bounded operator from 
H 2 (R;H 3 / 2 (doj)) into H 2 {M.;H 1 l 2 {duj)). 

3. Determining the twisting function from the DN map: an open problem 

Open problem. We first detail the reason why the problem of the identification of 8 from Ag remains 
open, at least in its full generality. To this purpose we consider a nonempty open subset 7 of duj and 
T = 7 x (-2L, 2L), for some fixed L > 0. Next we introduce the two following functional spaces 

H^H^iduj)) = {/ G i/ : (R; i/ 1 / 2 (9w)); supp/ C T}, 

and 

W(R) = <f(l); 0' G VF 1 - 00 ^)}!! 
Further, for i = 1,2, we pick ft, G W(R) and set A* = A(.t', ^(x 3 )), As = A 0j . In light of (|2.15|) , we see that 

(3.1) ((Ai - A 2 )ui,u 2 ) = / (Ai - A 2 )V«i • Vw 2 da;, 

for any function Ui G i? 1 (f2), i = 1,2, which is a weak solution to the equation div(AjViti) = in fl. Assume 
that 

01(33) = 02(33), |x 3 | > L. 

Then we deduce from (|3.ip that 

(3.2) ((Ai - A 2 )mi,m 2 ) = / {A 1 - A 2 )Vu 1 -Vu 2 dx 1 
where Q L = u) x (— L, L). 

Next, for < p < po, where po is some characteristic constant defined in jAG2| which depends only on ui 
and L, we put 

T p = {xe T; dist(a;,9r) > p}, 

U p = {x G R 3 ; dist(x,r p ) < p/4}. 

Upon eventually shortening po, we may assume without loss of generality that To = 70 x [— L, L] C T p for 
some 70 (s 7. Then, in view of |AG2j . we can find a Lipschitz domain Q p satisfying simultaneously: 

Q, C ftp, Tq C dfl n rip (e r and dist(x, Sfi p ) > p/2 for all x G U p . 

Moreover we know from jAGlj [Section 3] that there exists a unitary C°° vector field v, defined in some 
suitable neighborhood of du x (— 2L, 2L), which is non tangential to dfl and points to the exterior of il. For 
x° G r p , the point z T = x a + tv obeys Ct < dist(z T , 3f2) < r for all < r < To, according to jAGlj [Lemma 
3.3], where C and To are two positive constants depending only on ft, A, t and t. 

In light of [HKj . the operator div(AiV-), i = 1,2, has a Dirichlet Green function d = Gi(x,y) on ft p . 
More specifically Gi(x) = Gi(x,z T ), i = 1,2, is the solution to the boundary value problem: 

div(A,VG t ) = -8{x - z T ) in @'(fl p ), 
Gj = on dQ, p . 

Moreover, the claim of |AG2] [Corollary 3.4] remains essentially unchanged for the unbounded domain ft p 
arising in this framework: 



(3.3) 



Lemma 3.1. There exist two constants tq = To(fl,p) > and C = C(ft, \,t,t) such that the restriction 
Gi(.,z T ) to ft, belongs to H 1 (ft), and satisfies 

(3.4) \\Gi(;z T )\\ HHn) <Ct~ x / 2 , 0<t<t . 



2 We could choose W(R) = {9 G £>'(R); 9' G W 1 ' 00 ^)} as well. 
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Proof. Choose t = ro(0, p) so small relative to p that dist(z r , dil p ) > r for all r S (0, r ]. Since the ball 
centered at z T with radius t/3, noted _B r / 3 (z r ), is embedded in Q p \ O, we have 

(3.5) ||G i (-,^)llr=. 1 (n) < ||G i (-,« r )||y a ,i(n p \B T/ ,(« r )) < Cr" 1 / 2 , < r < r , 

directly from [HK : [Formula (4.45)]. Here C is some positive constant depending only on fl, A, t and i, and 
the space 

F 2 ' 1 ^) = {u e L 6 (n); Vu 6 £ 2 (ft) 3 }, 

is endowed with the norm ||u||y2,i (0) = ||u|| L 6( ) + || Vu|| L 2 (0)3 . 

Finally, Q having the Poincare inequality property since u> is bounded, we may actually deduce (|3.4[) from 
(US]). □ 

Now, as a Green function is a Levi function, it behaves locally like a parametrix. Hence, by applying 
[M] [Formula (8.4)], Gi, i = 1, 2, can be brought into the form 

(3.6) Gi(a;) - C(det (^(z,)))" 1 / 2 (A(z T ) _1 (^ - «r) ■ (a: - ^)) _1/2 + fl,-(a:), 
where C > is a constant and the reminder i?i obeys the condition: 

3(r ,a) et; x (0,1), ^(x)! + |z-z T | |Vi? 4 | < C\x-z T \- 1+a , x e Sl p ,\x - z T \ < r . 

Since 

div(A l VG 1 ) = in Q, 
in the weak sense, we can apply (|3.2[) ; whence 

<(A 1 -A 2 )Gi,G 2 )= / (A 1 -A 2 )VG 1 -VG 2 dx. 
Moreover, taking into account that Gi\ gn S i?p(K; H 1 ^ 2 (duj)), we obtain that 



/ (A 1 ~A 2 )VG 1 -VG 2 dx< - A^\\\\G 1 \\ m{R . m/ 2 {dLu)) \\G 2 \\ m{R . m 
Jn L 



>(9u,))) 



where Af, i = 1,2, is the restriction of Aj to the closed subspace H^(R; H 1 / 2 {dui)). As a consequence we 
have 



(3.7) / (A, -A 2 )VG 1 -VG 2 dx<C\\A T 1 -Al\\ \\G, \\m m \\G 2 \\ m 

Now, pick x° € To such that |6^(xg) — 9' 2 (x®)\ = \\9[ — #2 1 1 £°° (-£,£)■ Actually, we may assume without loss 
of generality that we have \9[{x%) - 9' 2 {x%)\ = 9 1 (x%) - 9 2 {x%). In view of (0[6]) and [XG2] [Formula (4.2)], 
the main term in the left hand side of (|3.7j) has the following expression 



[A-^jx )'^) - A-\(X )\S Q )] (X - Z T ) ■ (X - Z T ) 

3/2 3/2 ' 

B(z T , P )nn [Pq(x - z T ) ■ (x - z T )\ [Qo(x - z T ) ■ (x - z T )\ 

with to = 9' 1 {x%), So = 9 2 {x®), Po = A~ l (z' T ,to) and Qo = A^ 1 (z' T , sq), z t being the same as in above, i.e. 
z T = x° + tv for some r € . The main ingredient in the analysis developped in |AG2] is the ellipticity 
condition AG2J [Formula (2.5)] imposed on dtA(x',t). Indeed, this assumption entails 

,„ a , f [A-H(x°y, to) Sq)] (x - z r ) ■ (x - Z T ) 

(3.8) / 3^ JJJ dx - Grifo - so), 

JB(z T , P )nn [P (x - z T ) ■ (x - z T )] [Qo{x - z T ) ■ (x - Z T )\ 

for some constant C > 0, which leads ultimately to the desired result. Unfortunately, in this framework, the 
ellipticity condition fA~G2] [Formula (2.5)] is not fulfilled by d t A(x',t). This can be seen from the following 
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explicit expression 

Ai =0, 



A 2 = \x'\H- ^\x'\H 2 + \x'\ 



2 



A 3 = \x\ 2 t + ^\x'\H 2 + \x'\ 2 , 

of the eigenvalues of dtA(x', t), showing that the spectrum of dtA{x' , t) has a negative component for x' G du>. 
Moreover, it can be noticed that, due to the occurence of this negative eigenvalue, the weak monotonicity 
assumption |AG1| [Formula (5.7)] is not satisfied by the conductivity matrix under consideration either. 
Therefore, the approach developped in |AG2j does not apply to the problem under study. This explains why 
the determination of the twisting function from the corresponding DN map remains an open problem in the 
general case. 

Nevertheless, we shall establish in the coming section that this is not the case for affinc twisting functions 
anymore. Prior to examining this very peculiar framework, we first address the case of twisting functions 
which are close to some a priori fixed constant value. 
The case of twisting functions close to a constant value. Put 

(l + xi. —XiX\ -x 2 \ 
-x 2 xi l + x\ xi , x' G uj, t G M, 
-x 2 Xi I J 

and denote by A* g the DN map Ag where A*(x' ,9(x 3 )) is substituted for A(x' ,9(x 3 )). Then, by arguing as 
in the derivation of (|2.16[) . we obtain that 

(3.9) <(Afl-A5)u,u*)= [ (A(x',e(x 3 ))-A*(x',e(x 3 ))Vu-Vu*dx, 

for all solutions u,u* G iJ 1 (SI) to the equations 

div(A(x',9(x 3 ))\7u) =0 infi, 
div(A*(x',9(x 3 ))\7u*) = in Q, 

in the weak sense. Let us now assume that \\9' — l||oo < ^0- In light of ()3.9[) . we get that 

W^e - ^-e\\c(H- 1 {W,H- 1 / 2 (doj)),H 1 (R,H 1 / 2 (duj))) < C\\9' - lHoo, 

where C = C(w) is some positive constant. Therefore, if 9' is sufficiently close to 1, the operator A* g may be 
seen as a suitable approximation of Kg in the inverse problem of determining 9' from Ag. The main benefit 
of dealing with A g instead of Ag, boils down to the fact that 

d t A*(x',t) = A*(x', 1) = A(x', 1), for all and t G R, 

so the cllipticity condition AG2 [Formula (2.5)] required by the method developped in |AG2| . is fulfilled by 
A* . Therefore, by repeating the arguments of the proof of |AG2j [Theorem 2.2], we obtain the: 

Theorem 3.1. Let L > and M > 0. Assume that 6>i, 9 2 G W(R) obey 9' 1 (x 3 ) = 9 2 (x 3 ) for \x 3 \ > L and 
that 

|| 9[ || wn,°°(R), ||^2ll VK 2 .°°(R) < 

Then there exists a constant C = 0(M,u>, L) > such that we have 

\Wl - #2lU°°(R) < C||(^-ei) r _ (^8 2 ) T \\c(H- 1 (R,H- 1 / 2 (duj)),H 1 (R,H 1 / 2 {du)))- 



Note that 1 can be replaced by any constant. 
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4. The case of affine twisting functions 

In this section we address the case of affine twisting functions by means of the partial Fourier transform 
T X3 with respect to the variable x 3 . This is suggested by the translational invariance of the system under 
consideration in the infinite direction x 3 , arising from the fact that the matrix A(x',9'(x 3 )) appearing in 
()1.3jl is independent of x 3 in this peculiar case. 

In the sequel, we note £ the Fourier variable associated to x 3 and we write w instead of F X3 w for every 
function w = w(x' 1X3): 

w(x',£) = (F X3 w)(x' ',£), x' £ w, ^1. 

The first step of the method is to re-express the system (|1.3j) in the Fourier plane x' e £ R}. 

Rewriting the boundary value problem in the Fourier domain. We start with the two following 
useful technical lemmas. 

Lemma 4.1. For every w £ H 1 ^) it holds true that d Xj w = d Xj w, j = 1,2. 
Proof. Fix j = 1, 2. For every <p £ C^(uj) and -0 £ <S(R), we have 



<p(x') i / d x .w(x',x 3 )ii(x 3 )dx 3 ) dx' = / tjj(x 3 ) I / d x .w(x' ,x 3 )ip(x')dx' 

1 \JR / JR \Jw 

from Fubini's theorem. By integrating by parts in the last integral, we obtain 

d Xj w(x' ,x 3 )ip(x')dx' = — w(x' ,x 3 )d Xj (p(x')dx' , a.e. x 3 £ R, 



dX3 



so we get 

/ tp(x') I / d Xj w(x',x 3 )tp(x 3 )dx 3 ) dx' = - i>(x 3 ) I / w(x' ,x 3 )d Xj ip(x')dx' J cfa;3 

- -/a,, vM (f»(*>3)feWW. 

Further, the operator J- X3 being sclfadjoint in L 2 (R), it holds true that 

/ w(x' ,x 3 )-tp(x 3 )dx 3 = / uj(x',£)ip(g)d£, a.e.x'euj, 
Jr Jr 

whence 

tp(x') ( / d Xj w(x' ,x 3 )%i>(x 3 )dx 3 )dx' = - / d Xj tp(x') ( / w^', QipiQdt; ) dx' 
\Jr J Ju \Jr J 



by integrating by parts. From the density of C^°(o;) in L 2 (oj), the above identity entails that 

d Xj w(x' ,x 3 )ip(x 3 )dx 3 = / d Xj €i(x' ) ^)ilj(^)d4, a.e. x' £ cj, 



for every j/> £ <S(R). From this, the selfadjointness of and the density of <S(R) in L 2 (R), then follows 
that 9^(0 = d Xj w. □ 

Lemma 4.2. Let C = (Cki)i<k,i< 3 £ W 1 ' 00 (oj) 3x3 be such that (Cm {x'))i<k,i< 3 * s symmetric for any x' £ lj. 
Then every w £ iJ 1 (SI) obeying 

(4.1) [ CVw-\7vdx = for allv eH^fl), 
satisfies the equation 

(4.2) - div x , (C^OV^u;) +P(x',£) -V x >w + q(x', £)w = mZ?'(ft), 
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where 



C(x') = (Cy(l'))l<^<3 



q(x',0 = -i£div x , +ec 33 (x>). 

Moreover, if w G H 2 (Q) is solution to (|4.1j) £/ien £/ie identity ()4. 2|) ZioWs true for a.e. (x',£,) G fi. 
Proof. Choose u = y> ® in ()4.1[) . where y> G Cq°(uj) and ?/> G so we have: 

(4.3) I C kl (x')d Xk w(x',x 3 )d Xl {ip®$)(x',x 3 )dx'dx 3 = 0. 

fe,/=l,2,3 n 

For 1 < fc, Z < 2, we notice that 



C k i(x')d Xl w(x',x 3 )d Xk ((p<gi : tp)(x',x 3 )dx'dx 3 = / C fe £(x')9 a;(i; ^(x') I / d Xl w(x' ,x 3 )i>{x 3 )dx 3 ) dx' 

Ju> \Jm. J 

(4.4) = f C H {x')d xl w(x',ti)d Xk {y®i))(x',Odx'd(i, 

Jn 

directly from Lemma I4TT1 Further, for all I = 1,2, it holds true that 

C 3 i(x')d Xl w(x',x 3 )d X3 ((p<^i))(x',x 3 )dx'dx 3 = / C 3j (x')(p(x') I / d x ,w(x', x 3 )$'(x 3 )dx 3 ) dx'. 



C 3 z(x')y>(x') ( / 9 2;i ti;(x',X3)(-^)V'(x3)dx3 ) dx'. 
Therefore we have 

(4.5) / C 3 i{x l )d Xl w{x',x 3 )d X3 {ip®^){x , ,x 3 )dx'dx 3 = / C 3 i(x')(-i£)d Xl w(x' , ® 4>)(x' , £)dx'd£. 
Jq Jn 

Next, for k = 1,2, we may write that 

C k3 {x')d X3 w{x' ,x 3 )d Xk (ip® tp)(x',x 3 )dx'dx 3 = / C k3 {x')d Xk ip(x') / d X3 w(x' , x 3 )ii{x 3 )dx 3 I dx' 
! Jut \Js, J 

= C k3 (x')d Xk ip{x') ( / w(x',x 3 )$'(x 3 )dx 3 ) dx', 

Jul \JR J 

with 

/ w(x' , x 3 )$ (x 3 )dx 3 = - / w(x',x 3 )(i€)ijj(x 3 )dx 3 = - / (i£)w(x', £)il>(£)d€, 
it ii Jk 

hence 

/ C k3 {x')d X3 w{x' ,x 3 )d Xk {ip®Tp){x' ,x 3 )dx'dx 3 = I C k3 (x')(i£)w(x' , £)d Xk {if ® ^)(x', £)dx'd£. 
Jn Jn 

By integrating by parts in the right hand side of the above identity we obtain that 

C k3 {x')d X3 w(x',x 3 )d Xk ((p®'$)(x',x 3 )dx'dx 3 = \ C k3 (x')(-i£)d Xk w(x',£)((p <g> ip)(x',£)dx , d£ 



(4.6) + / d Xk C k3 (x')(-i0w(x',Z)(cp®iP)(x',0dx'dt. 

Jn 

Finally, bearing in mind that d X3 ((p (g> if)) = ip ® (—i£)ip and noticing that 

4.»(x',, 3 ,(^(, 3)fc . - 1 » (I ', I3 )(^fe)^ - - ^.ohomo* 
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we find out that 

(4.7) / C 33 (x')d X3 w{x\x 3 )d xs {ip®i>)(x',x 3 )dx'dx 3 = / C 33 (x')(-i£) 2 w(x', £)(<p <g> ff>)(xf, 

Now, putting (|4.3p - (|4.7p together, we end up getting that 

(-div x ,(C(x')V x ,Q) + P{x',0 ■ V x ,w + q(x',Ow,$) = 0, $ G C °°H ® C °°(K), 

where (•, •) denotes the duality pairing between C£°(f2) and P'(fi). From this and the density of Cg°(ui) ® 
C£°(R) in C£°(ft) then follows that 

-div x /(C(a/)V a /«;) + P(a/, • V x ,w + q(x', t)w = in V'{Q), 
which completes the proof. □ 

In the remaining of this section we assume that #(£3) = ax 3 + 6, where a and b are two fixed real numbers. 
With the help of Lemma \A. 2 1 we may now re-express (|1 .3[) in the Fourier plane. For the sake of simplicity, 
we shall write A a (x') instead of A(x', 9(x 3 )), as we have 9'(x 3 ) = a for all £3 G ML 

For g G ff x (M) n L 1 (K) such that J R g(x 3 )dx 3 = 1 and for h G # 1/2 (&j), we consider the /f 1 (Q)-solution 
u to (|1.3|) . with: 

f(x',x 3 ) = g{x 3 )h(x'), x' G duj, x 3 G R. 
Since u is solution to P~Tj) with C = A Q G W l >°° (u;) 3 * 3 , we deduce from Lemma |4~21 that u G L 2 (R; if 1 ^)) 
is solution to the system 

-div x ' (A a (x')V X 'u(x', ()) - 2ia£x ,A - • W x >u + ( 2 u = in £>'(0), 

() = £(()/ on Ski;, for all (el, 



(4.8) 
where 



x = (—X2,xi) and A a (x') 



1 + x 2 a 2 —xix\a 2 
—X2X\a 2 1 + x\ a 2 
We turn now to examining the variational problem associated to (|4.8[t . 
Variational problem associated to (|4.8p . Let us consider the bilinear form 



A i [(v,w),((p,ip)]= / A a Vv ■ Vipdx' - 2a( / x'- 1 • Vwipdx' + ( 2 / v<pdx' 

+ / i„Vw • Vipdx' + 2a£ / x'- 1 - • Vv^dx' + ( 2 / wtpdx', (v,w),((p,ip) eH, 



defined on the Hilbert space H = Hq(lu) x Hq(u>) endowed with the norm 

||( W , W )||„ = (||V U || 2 2M + ||V W || 2 2M ) 1/2 . 

Taking into account that 

A a {x')C, ■ C > |C| 2 , for all C G M 2 and x' G u, 

and that 

2a|£| 





x ■ Vvw 


dx' < a 2 S 2 f 






J UJ 



where 5 = m&x x i etJ \x'\, it is easy to see that 

(4.9) Ad(v,w),(v,w)}>(l-a 2 S 2 )\\(v,w)\\ 2 n . 

Let us fix ao > so small that a = 1 — Oq<5 2 > 0. Then, due to the above estimate, the bilinear form A^ is 
a-elliptic for every (el, provided we have \a\ < ao. For each $ G C(R; H') and every (el, there is thus 
a unique (v(£),w(£)) G H satisfying 

(4.10) MW0M0)>(<P,1>)] = (*(0,(<P,1>)) forall (<p,tl>) eU, 
by Lax-Milgram's lemma. From this then follows that 

A+vMZ + v) - w(0, + v)- M0)> 0= V 1 )] 

(4.11) = ^ [(«(£), «;(()), (p, V^)] - ^W(»(0, VO] + <*« + »7) - *(0. (<P, ^)). 
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for each £, r/ G R and (95, ?/>) € Further, by noticing through elementary computations that 
A^ +n [(v,w),((p,ip)] = A^ +ri [(v,w),((p,rp)]-2ari / x' 1 - ■ (y>V x ' to -ipV x <v) dx' + 77(2^ + 7?) / (vip + wip)dx' , 



for every (<p,ip) G H, we deduce from (|4.11[) and Poincare's inequality that there is a constant C = 

C(£,uj 7 a ) > satisfying 

■4* [MOXO), M£ + »?) - «(0. + »?) - «(0)] 

- Az +n [(v(0,w(0), + »7) - «(0> «>(£ + - «>(0)] 

(4.12) < cMUKO.^IkllK^ + ^-^O.^ + ^-MOJIIw. 

for all £ 6 R and rj G [-1, 1], In light of ([3~5]t. P~TT|) written with (<p,ip) = (v(£ + rj) -v(£),w(£ + r)) 
and (|4.12p , wc thus find out that 

This proves that (v,w) G C(R;H). Moreover, we obtain 

(4.13) IIMO>w(Olk<(l/a)ll*(OII*', 

directly from (j43| - (|^TU|) . Further, it is easy to check for $ G C^R;^') that (</(£), w'(£)) 6 C(R,U) is the 
solution to the variational problem 

^[((«'(0X(0), foVO] = (*o(0. falM + ($'(0, (V,^)) for all ( V> V) G H, 
where we have set 

<*o(0.(^)} =2a / ^ • (pV^Cf) - VV^(0)&/ + 2£ f (v(£)cp + w(t)iP)dx'. 

J to J to 

Using (|4.13p and putting {£) = (1 + l^l 2 ) 1 / 2 , we deduce from the above estimate that 

\\(v'(o,iv'(0)\\n < c((o\\m\\w + \mo\\w), £ e R, 

for some constant C = C(ao, w) > 0. Similarly, if $ G C 2 (R; W) then the same reasoning shows that(u, «;) G 
C 2 (R,-H) and 

ik«"(o,«;"(o)IIk < c«e) 2 ii*(oiiw + <oii*'(oii«' + wm\n>), c e ». 

Hence we have obtained the: 

Proposition 4.1. For every $ G if 2 swc/i that (0 2 ~ J $ (j) G L 2 (R,H')- J = 0, 1, the variational 

problem (|4. 1 0[) admits a unique solution (v,w) G H 2 (M.;H) satisfying 

2 

(4-14) \\(v,w)\\ H 2 {Rm < C(J2 Utf-^WisW)), 

0=0 

for some constant C = C(uj,ao) > 0. The above assumptions on $ are actually satisfied whenever $ = ^ for 
some * G H 2 (M.;W) such that x 3 $> G H 1 and x^ G L 2 (R;H')- Moreover, the estimate gHH) reads 

2 
0=0 

in this case. 

Armed with Proposition 14. II we may now tackle the analysis of (|4.8[) . 
Analysis of the solution to (|4.8[) . Pick F G H 1 ^) such that F = f on duj and H-FH^-i^) < C(cj)||/|| J3 -i/2( aw -). 
Let u r (resp. u l ) denote the real (resp. imaginary) part of u = u — = + it?. Since the Fourier 

transform u of the i7 1 (il)-solution u to p. 31) is actually solution to (|4.8p . we get by direct calculation that 
(u r ,u l ) is solution to the variational problem (|4.10[) . where 

(4.16) m),M)) = -M(S r F,g i nM)], 
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and g r (resp. g 1 ) stands for the real (resp. imaginary) part of 5. Further, in light of (|4.16|) we check out 
using elementary computations that 

(4-17) WHOWw < C(0 2 \g(mf\\m^) 

(4.18) < C((0 2 \g , m + (0\m\)\\\f\\Hy H a.) 

(4.19) wm\w < c((o 2 \rm + W(o\+m)\)\\f\\ H ^a u) , 

for some constant C = C(a ,w) > 0. Therefore we have (£) J '$( 2-J ') G L 2 (M) for j = 0,1,2, provided 
(0 y 3 "= -k 2 (^)i this later condition being ensured by the assumption x 3 3 g G H 4 ~ 3 (M.). From this and 
Proposition 14 . 1 1 then follows the: 

Corollary 4.1. Assume that g G H 4 (R) is such that x 3 g G H 3 (R), x\g G H 2 (R) and J R g(x 3 )dx 3 = 1. Then 
it holds true that u G H 2 (R; H 1 (a;)) . Moreover we have u € L 1 (R; H 1 (u>)) andU = u(-,0)= L 2 3 )da;3 G 
if 1 (w) is the variational solution to the following boundary value problem 

{div x > (AaVa't/J = inus 
U = f on dtJ. 

In view of (j4~5|) and (gTTUJ), we deduce from (t4~TT)) - (l4~Tm) that 

||di Va: , (A a V^=u(-,0)llx-(^) < CiO^giOl WfWm^au.) 
\\dsdiv x ,{A a V x >u(;Z))h*M < C((0 5 \g'm + (0 4 |3(0l)lll/ll^/ 2( a.) 

ii9|div^(i a v K m(,o)iu 2M < c({o 6 m)\ + (o 5 \d'm + (o 4 i5"(oi)ii/n ffl / 2(aw) , 

for some positive constant C depending only on ao and uj. From this and Corollary 14.11 follows the: 

Proposition 4.2. Let g G H 6 (R) be such that x 3 g G H 5 (R), x 3 g G H 4 (R) and J R g(x 3 )dx 3 = 1. Then we 
have A a V x ,u- v(x') G H 2 (R; H~ 1 / 2 {du>)) and thus A a Vu-v{x) G L X (R; H~ 1 l 2 {du:)), with 

A a V x ,U ■ v{x') = A a V x >u{; 0) • v{x') = [ A a Vu(; x 3 ) ■ v{x')dx 3 G H~ 1 ^ 2 (8uj). 

Jr 

In light of Proposition 14.21 we now introduce the two following DN maps: 

A a : / G H 1/2 {du) A a \7u{-, x 3 ) ■ u(x) G L\R; H- 1/2 {du)) 
A a : / G HV*(du) -> A a V x ,U ■ v(x') G H- 1 ' 2 {duj). 
These two operators are bounded, and they satisfy the estimate 

(4-20) l|Ai -A 2 || £(H i/ 2( a w ),H-v=(9a;)) < II Ai - A 2 \\ c ( H1/a{du)tLl0t . H - wa(du ))) ' 

where, for the sake of simplicity, we write Aj (resp. Aj) for A aj (resp. A a .), j = 1,2. 

Finally, the last step in the analysis of the case of affine twisting functions involves noticing by direct 
calculation that the matrix d a A(x',a) has two eigenvalues Ao = and Ai = \x'\ 2 . Hence, by mimmicking 
the proof of [XGl] [CLAIM, page 169], we obtain the: 

Theorem 4.1. Let ao > and let g be the same as in Proposition \4-.2\ Assume moreover that 1 — a^S 2 > 0. 
Then there exists a constant C > 0, depending only on ao and 10, such that the following stability estimate 

\ ai -a 2 \ <c\\^i-^4^ m/Hdu)tL1{u . H - 1/Hdu)) y 

holds true whenever 

\a>i\, |a 2 | < a . 
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5. The DN map for the original problem 
We first start by defining the trace space for functions in iJ 1 (f2e). To this purpose we set for all L > 0, 

fig = {(Re( X3 )x',x3); x = (£1,2-2) e u, x 3 e (-£, L)} 

and 

r e = {(Re(x 3 )X , x 3 ); x = (xi,x 2 ) G <9w, £3 G [— L, L]}. 
Given L > 0, it holds true for every u G H 1 ^) that u G H 1 (VLq), whence u^l G H 1 / 2 ^^). Putting 

# 1/2 (r£) = {h = g [T L in L\T L e )- g G i/ 1/2 (^)}, 

we thus have U\ t l G i7 1 / 2 (Fg). Here H x / 2 {Tg ) is equipped with its natural (quotient) norm 

ll /l llHi/2 ( r^) =inf{||s||jjV2(an£); 5|r£ = M- 

Further we define 

H^(dQ e ) = {h G Lf oc (^); ft|rf 6 # 1/2 (r£) for all L > 0}, 
1/2 

and introduce the following subspacc of _ff lo ' c (<9f2g): 

H 1/2 (dn ) = {he H^(dCl e ); there exists v G i? 1 ^) such that v\ 9Ug = h}. 
Here and henceforth v\qq b = h means that v^t — h^L in the trace sense for all L > 0. It can be checked 
that iJ 1 / 2 (9ile) is a Banach space for the quotient norm 

IMIsva(en e ) = inf {IMIi? I (n«); w |an e = M- 
Let us now introduce the mapping 

.9 — > f = g°fe, 

where, for the sake of shortness, we note 

ipe(x) = T e{x3) (x',x 3 ). 

Further, pick g in C^dflg) and let v G Cq(IR 3 ) be such that v\gn g = g. Setting u = u|q 8 oyjg, we notice that 

\\Ie9\\m(R;W/2(du,)) <C{uj)\\u\\ H i {n) < C(u, 6)\\v\\ H i( Qe y 

Therefore, we have 

(5-1) \\hg\\m(R;H^Hdu)) < C(u>,0)||$r|| 5l/a(ane) for an y 9 G Co( 5n «)- 

Let us now consider g G if 1 / 2 (9r2e) and w G H 1 (ttg) such that Wian e = ,9. For any sequence (v n ) n G C*d(R 3 ) 
such that v n \Q 9 — > v in H 1 (flg), it is clear that 

lis - 5n||jfi/2( ane ) < ll« - «T»||Hi(n e ). 

provided we have g n = i>n|an 9 - Hence (g n )n converges to g in H l / 2 {d£le) 

For all n > 1, put /„ = 7 p„ = o <^ e and u„ = w„ o <^ e . Since /„ = u n \ dn , we see that 

||/n - /m||ifi(R;ffi/2( 9tJ )) < C(w)||w n - tin, II ffi(n) < C(n,0)\\v n - U m ||jyi(n e ). 

Consequently, (/„)„ is a Cauchy sequence in iJ x (R; H l / 2 (du>y) so there exists / = lim„ /„ G H 1 (R; iJ 1 / 2 (9w)). 
Set / = leg. In view of (|5.1[) , Ig extends to a bounded operator, still denoted by Ig, from H 1 ^ 2 (dfte) into 
iT^R; if 1 / 2 (9a;)). 

Arguing as above, we find out that the mapping 

Jg : Cl{dti) ^ CUdttg) 

f — > 9 = f ° i>e, 
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where tpg = (fig 1 , extends to a bounded operator, which is still called Jg, from iJ 1 (R; H 1 / 2 (du>)) into 

H^Hdng). 

Evidently, we have IgJgf = / for all / G Cq(90) and Jglgg = g for all g G C^dflg). Therefore Jg = Ig l , 
by density. 

Next, by reasoning in the same way as in the derivation of (|1.3|) . we prove with the help of the Lax- 
Milgram's lemma that the boundary value problem has a unique solution v G H 1 (ilg) for every g G 

H^idfle). Mor eover the operator Ag is well defined as a bounded operator from H^ 2 (dn e ) into its dual 
space H-^idfLe). Similarly to Ag, it is can be checked that Ag is characterized by the identity 

(Agg, h) = I Vu ■ VHdy, 

which holds true for all h G H x ^ 2 {dVlg) and H G if 1 (fig) such that H\qq b = h. By performing the change of 
variable y = fg(x) in the last integral, we get that 

(A e g,h)= / AVu- V(H oifg)dx, 
Jn 

where u is the solution of the boundary value problem (jl.3p with / = Igg. Therefore we have 

(A e g,h) = (Aglgg,lgh), 
which means that Ag = IgAglg, or equivalcntly that Ag = JgAgJg. 
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